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Abstract 

In this paper we formulate a new iV-state spin integrable model on a 
three-dimensional lattice with spins interacting round each elementary cube 
of the lattice. This model can be also reformulated as a vertex type model. 
Weight functions of the model satisfy tetrahedron equations. 
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1. Introduction 



Recently in Ref. [1] Korepanov has constructed a new solution of tetrahedron 
equations with spin variables lying on the edges of a 3D cubic lattice [2]. It 
leads to a commuting family of transfer-matrices and possible integrability 
of this vertex model on the cubic lattice. Corresponding spin variables take 
N — 2 values. There are only 16 nonzero weights Rjli^i^ from 64 possible 
ones. 

Later Hietarinta in Ref. [3] has proposed another vertex solution with 
16 nonzero weights. It strongly reminds a solution of Ref. [1], but four 
weights from sixteen in these two models are different. This model satisfy to 
a duality B property in a terminology of Ref. [4] and can be reformulated as 
an interaction-round-cube (IRC) model with weight function W{a\efg\bcd\h) 
depending on eight surrounding spin variables of an elementary cube of the 
lattice. Note that Baxter in his paper [5] also has reformulated Zamolod- 
chikov model [6] with spin variables belonging to faces of 3D lattice as an 
IRC model. We fail to generalize the original Korepanov solution for an 
arbitrary N. 

In this paper we propose a generalization of Hietarinta solution of tetra- 
hedron equations for an arbitrary number of spin variables. The obtained 
solution has a simple multiplicative form, depends on four parameters, which 
have a simple geometric interpretation. Namely, one can choose as param- 
eters of our solution four angles of an arbitrary quadrilateral with two di- 
agonals. Hietarinta solution appears in the case if this quadrilateral can be 
inscribed in the circle and N = 2. Also note that in the case N = 2 our 
solution can be obtained as a particular case of Zamolodchikov model [5,6]. 

The paper is organized as follows. In section 2 we recall a formulation of 
IRC type models , introduce necessary definitions and propose an ansatz for 
weight functions. In section 3 we show that the tetrahedron equations are 
reduced for this case to some special identity, which is proved in appendix. 
Section 4 is devoted to a parameterization of obtained A^-state solution. And, 
at last, in section 5 we discuss briefly results of this paper and its possible 
generalizations. 
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2. Formulation of The Model 

In this section we recall some definitions following Refs. [7,8]. Consider a 

simple cubic lattice L and at each site of L place a spin variable taking its 
values in Zat, for any integer N >2 (elements of are given by N distinct 
numbers 0, 1, . . . , A*"— 1 considered modulo N). Allow all possible interactions 
of the spins within each elementary cube. The partition function reads 

^ = E n W{a\eJ,g\h,c,d\h), (2.1) 

spins cubes 

where W {a\e, f , g\b, d\h) is the Boltzmann weight of the spin configuration 
a, . . . , /i. 

We need in some notations to formulate our ansatz for the Boltzmann 
weights. Denote 

uj = exp(27ri/A^), o;^/^ = exp(7ri/A^). (2.2) 

Taking a;, y, z to be complex parameters constrained by the Fermat equa- 
tion 

+ = (2.3) 

and I to be an element of Z^, define 

-fej/.-IO-riTZ^- (2.4) 

This function has a following property 

w{x, y, z\l)w{z, uj^^'^y, ujx\ - = 1, I e Z^, (2.5) 

where 

$(0 = (2.6) 
The tetrahedron equations can be written as (Eqs. (2.2) of Ref. [5]) 

^ W{a4\c2, ci, C3I61, 63, b2\d)W'{ci\b2, 03, 6i|c4, d, celb^) 

d 

X W"{bi\d, C4, C3\a2, 63, b4\c5)W"'{d\b2, 64, &3IC5, C2, cejai) 
= E C4, C3|a2, 04, a3\d)W"{ci\b2, 03, a^ld, C2, Cejai) 

X W'{a4\c2,d,C3\a2,b3,ai\c5)W{d\ai,a3,a2\c4^,C5,CQ\b4), (2.7) 
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where W, W, W" and W" are some four sets of Boltzmann weights. As has 
been shown in Refs. [2,9], this relation ensures commutativity of layer-to- 
layer transfer matrices constructed from W and W weights: 

T{W)T{W') = T{W')T{W). (2.8) 

We are looking for a solution for (2.7) in the following form: 

W{a\efg\hcd\h) = uj^h-f)(a-h-e+h)^ 



w(p2,Pi2,Pi|a -h-e + h)w{p4,ps4,ps\ 


\— a + c + f — h) 




\-b + c-e + f) 



(2.9) 



where functions w are defined by (2.4) and coordinates Pi, Pij satisfy Format 
constraint (2.3). We imply that the set of weight functions W depend on 
parameters p^, p^ - and etc. 



3. Proof of the Tetrahedron Equations 

In this section we will show that tetrahedron equations (2.7) for weight func- 
tions (2.9) are equivalent to the identity (A. 8) provided that parameters of 
weights satisfy some special algebraic constraints. 
Substituting formula (2.9) in (2.7) we obtain 





h-a) w{pI,pI^,pI 


|c'-6') 


1 III III III 


\c- d)w{pQ,pr,Q,p^\ 


d' — a') 



w Hp6>P56>P5I - a - a)w{p'4,p'^4,p'^\ - h - a)w{p'!l ,p'^4,p'I^\c + a) 

X«^(P2">Pl2>Pl'l 0-)«'(P4,P34,P3| - fc' + 'y)w{p2,p'i2^p'i\ - d' + a) 

xw{p2,Pi2:Pi\b' + d' -a' - (t)w-^{p'^,p'J^(^,p'^\c' - b' - d' + (7)00"^ ^^"-^'^ 



w{p'ip'U,p'i 


\c-b) W{p2,p[2,p'l\ 


b' - a') 


W{P6, P56,P5\d~ a) w{p'^, p'^'q, 


\c'-d') 



X (3.1) 



^ Hp6>P56>P5I - a' - (T)w(p4,p^j4,p^| - 6' - (x) w (p^" , , p'3" | c' + a) 

X«^(P2'>Pl2>Pl'l -d' - (7)w(p4,P34,P3| -b + Cr)w{p'^, p'(2: Pi \ - d + a) 

xw{p2,pi2,Pi\b + d-a- a)w-\plp'U,pl\c -b-d + (7)a;^'+^(^'-^), 
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where we made a replacement a ^ —a in the LHS, cr — >■ Ci + C5 — cr in the 
RHS and introduced new spin variables: 



a — —as + 62 + C4, a' — a2 — bs — Cl + C2 — C5, 

b = -as + 64 + Cl, 6' = 04 - 63 - Cl, 

C = — Oi + 64 + C2, c' = 04 — foi — Cl + C4 — C5, 

d = -ai + 62 + C5, = a2 - &i - C5. 



(3.2) 



Relation (3.1) looks very similar to identity (A. 8) from appendix. But 

for complete coincidence of spin structure of (3.1) and (A. 8) we need to use 
(2.5). Using (2.5) several times we can easily reduce (3.1) to identity (A. 8) 
with the following identification: 



Since we can choose four vectors x, y, z and u in an arbitrary way, we come 
to the four-parametric solution of tetrahedron equations (2.7). In the next 
section we will show that this solution can be parameterized by the angles 
of a degenerate tetrahedron with all apices lying in one plane. 

4. Parameterization of solution 

Using formulas (A.5-A.6) from appendix A we can exclude parameters en- 
tering in X, y, z and u from relations (3.3). Then we come to the following 
constraints among parameters p: 




(3.3) 



{xiy3,C,i^X3yi) ~ {p[,UJ^/^p[2,UJP2) 
{zsU„lziUs)r^iptu^'^'p'L^p'^) 



{y3Zl,C",UJyiZ3) - (P4,P34,P3), 
{xsUi,i",XiUs) ~ (P6,P56,P5)- 



P6 _ P2P4 p'e _ P2P4 Pe _ P2P4 Pq 




(4.1) 



' / / / ' // n in II, 
P5 P1P3 P5 P1P3 P5 P1P3 Pb 



P'2 _ P2P2 P'i _ p'aPI' P4P2 _ P4P2 



(4.2) 



Pi PlPl P3 P3P3 P3P1 P3P1 
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1 



PlP56 P3P12 P2'Pu 

P5P12 P2PU Ps'Pu 



1 



P5P34P3P56 P6Pu P5 
PiP56 P'5Pu pIp'U P& 



1. (4.3) 



P',P\2P'2P'kPyi2 



Also we obtain additional constraints on parameters pi corresponding to for- 
mulas (A. 6) for ^'^ and ^'^ from appendix A. These constraints ensure a 
periodicity {modulo N) of all w functions entering in (2.9) and can be easily 
obtained as a consequence of (4.1-4.3). 

Remind that the function 'w(j)j,p,ij,pi\n) depends on one complex param- 
eter Pi/pj and parameter p^j defined by a relation pf- = pf —pf can contain 
an arbitrary multiplier o;"'^, where riij e Z. 

It appears that instead of variables Pi, Pij and pj it is more convenient to 
introduce "angle-like" variables ai, 02 and 03 as 



P2 = exp(— ifli/A^), pi = exp(mi/A^), pu = {2ism{ai)Y^^uj'^^'^, 

P4 = exp(-m2/iV), p3 = exp(m2/iV), P34 = {2i sin(a2)}^/^u;"^%(4.4) 

P6 = exp(-ia3/A^), p5 = exp(ia3/A^), = {2i sin(a3)}^/'^a;"56 



Then relations (4.1-4.2) are reduced to linear constraints among Oj: 

— ai + 02, Og = a[ + 03, 03 = a" + Og, a'^' — a"' + 03", 
ai - a[ + a'l = 0, - a'^ + aHl = 0, 02 - 4 - a!^ + a'/' = 0. (4.5) 

Now let us fix a choice of phase multipliers a;"'^ . Taking into account 
formula (2.9) it is easy to see that we have the following dependence of the 
weight function from Uif 



In fact, this multiplier corresponds to a simple gauge transformation of the 
weight function and hereafter we will omit all niTiltipliers uj'^^k 
Then relations (4.3) can be rewritten in the following form: 



and etc 



,ni2(a-&— e+/i)+n34(-a+c+/-/i)+n66 (6— c+e-/) 



(4.6) 



sin(ao) sin(a^/) sin(a3') 
sin(ai) sin(a3) sin(a']") 



1 



sin(a3) sin(a'/) sin(a2") 
sin(ai) sin(a2) sin(a']") 



1 



sm(a2j sm^Ogj sm(a2j 
sin(a3) sEi(a^ sin(af) 




(4.7) 
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Formulas (4.5,4.7) have a nice geometrical interpretation (see Fig. 1) 



pill 




Fig. 1 

Then formulas (4.5) are valid in an evident way and formulas (4.7) are the 
consequences of sine theorem for four triangles formed by two sides and one 
of the diagonals of the quadrilateral. 

5. Discussion 

In the previous section wc obtained a four-parametric solution of the tetra- 
hedron equations for which each weight function depends on two spectral 
parameters. 
Namely 

W = W{ai, 02), W = W{a[, 4), W" = W{a'l, 4), W" = W{a"', a'^). 

(5.1) 

Weight functions are defined by multiplicative ansatz (2.9) where parameters 
Pi and pij can be calculated through angle variables ai, 02 from formulas (4.4) 
(remind that 03 = oi -|- 02). All set of angle variables a^, a^, a'l and a'l' is 
described by Fig. 1. 

It is naturally to ask: is this solution of tetrahedron equations new or not? 
To clarify this situation we have considered the only known solution of tetra- 
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hedron equations with N > 2 [10] corresponding to the model proposed by 
Bazhanov and Baxter [7] in the "plane" limit for which all apices of tetrahe- 
dron belong to one plane. More precisely, let us consider a parameterization 
of Ref. [8] and take the limit Oi, 62, 9^, 9^ ^ and 6*3, 9q it. This limit 
corresponds explicitly to parameterization (5.1) and to Fig. 1. Formulas for 
weight functions look as 



Xi = exp(ia2/A'")(sinai)''-/-^, X4 — u ^/'^ex.p{—iai/N){sma2y^'^, 
X2^exp(iai/N){sma2y^^, X3 = u;^/^ exp(-ia2/iV)(sinai)V^, (5.3) 
xu = (sin(ai + 02))^/^, X23 = (sin(ai + 02))^/^, 



and angles Oj appear as limit values of plane tetrahedron angles. The sub- 
script "0" after the curly brackets implies that the expression in the curly 
brackets is divided by itself with all exterior spin variables equated to zero. 
If the signs of a in w functions entering in (5.2) were equal we would use for- 
mula (A14) of appendix A of Ref. [10] and come to a multiplicative solution. 
It can be done for the case N — 2 where all spin variables belong Z2 and we 
can change all signs in an appropriate way. Only in this case we obtain from 
(5.2) our multiplicative solution up to some gauge transformation. 

So it seems that multiplicative solution (2.9) should be new for > 3. 
Unfortunately, we did not succeed yet in a generahzation of this solution to 
a five-parametric one which could be described by the angles of the usual 
tetrahedron. 

Note also that solution (2.9) can be naturally rewritten in terms spin 
variables lying on the edges of the lattice. 

Let us define the following spin variables [4] 




(5.2) 



where 




(5.4) 



(5.5) 
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Note that spin variables ii, i2, is, ji, j2, js should satisfy two constraints 



(5.6) 



Then tetrahedron equations take the vertex form 



ETjkiMM Tj'jikiks r>"j2jike D^Oajsje 
11,12,13 feiMis k2k4i6 ksk^ke 



ErtiiiksMM TD'/kikije TD'kijijs rihhh 
»3,«6,«6 1214^6 nfe4fe6 fclfe2fc3 

fcl.fe2.fe3. 
fc4.fc6'fc6 



(5.7) 



At last let us suppose that the quadrilateral shown in Fig. 1 can be 
inscribed in the circle. In this case all relations constrained angles Oj, a^, 
a'l and a'l' appear to be linear, and we have only three independent angle 
variables. For N = 2 this model coincides with that of Ref. [3] (to be exact, 
we are to change else aj — >• iaj). 

Acknowledgements. We are indebted to J. Hietarinta for giving an 
opportunity to get acquainted with his solution prior to publication and for 
the valuable discussion. Yu.G. Stroganov is also grateful to the staff of the 
Physics Department of Oklahoma State University, where a part of this work 
has been done, for their kind hospitality. 

Appendix 

In this Appendix we will prove the identity (3.1), from which it follows the 
tetrahedron equation for the multiplicative ansatz. 
We begin with the following relation: 




w{xiy3,^,ux3yi\a - b) 
w{z3Ui,^, ziusld- c) 



(A.l) 



where 
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z! = {UiZ2,i,U2Zi) 

u' = {xiy2,^,uJX2yi), (A.2) 

the vectors x, y, z, u are arbitrary, and the auxihary parameters ^ and ^ obey 
the relations 

e = x^y^ - xM 

= < - zi'u^, (A.3) 

phases of ^ and ^ can be chosen in arbitrary way. The subscript "0" after 
the curly brackets implies that the expression in the curly brackets is divided 
by itself with all the exterior spin variables equated to zero. 

This relation is equivalent to the relation (B9) from the Appendix B of 
Ref. [10] and can be interpreted as some symmetry transformation of the 
Boltzmann weights of the generalization of the Bazhanov - Baxter model, 
considered in Ref. [11]. 

Applying this transformation triple, we obtain the following relation: 



E 



w{x\a + (7)w{z\c + a) 







^atzj, w{y\b + a)w{u\d + a) 

^^n{b+d)+d{a-b+c-d) w{x3y2ZsU2,^',X2y3Z2Us\-n) 

w{xiy2ZiU2, ujX2yiZ2Ui\a — h + c — d — n) 



w{xiy3,C,ujX3yi 


\a - b) w{y3Zi, u}yiZ3\c - b) 


w{z3Ui,i, ZiUs] 


d-c) W{X3UI,1" ,XiU3\ 


d — a) 




w{x"' 


1 - c + a)w{z"'\ 


— a + a) 


w{f'\ 


-d + a)w{u"' 


\-b + a) 



\ 



where 



^a{-a+b-c+d+n) I (A.4) 
) 



x"' = {z3i',eix2,z,e) 
f' = iu3e/co,eiy2,u,e) 

^" ={x3l'A"^Z2,X,i') 

u"'^{y3i',C"Cu2,u;y,C'), (A.5) 
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and 



S — -^2 Us ^2 "3 



•^3 i/2 ^3 "2 



S — -^2 Ul ^2 -^1 y2 ^1 "2 

s — yi ^3 y3 ^1 



(A.6) 



Note that in the special case when xiy^ziUs = x^yiz^uiuj the summands 
from both parts of (A. 4) become the Boltzmann weights for the Bazhanov - 
Baxter model [7], and (A.4) becomes the well known Star - Star relation. 

Consider now the second copy of (A.4) with the same x, y, z, u, with the 
same number n and with some spines a',b',c', d' obeying 

a - 6 + c - d = a' - 6' + c' - d'. (A.7) 

Multiplying left hand side of one copy of (A.4) by the right hand side of the 
other copy, we obtain the identity, from which the w - functions containing 
the spin n are cancel out. Summing the obtained relation by n, we obtain 




w{x\a + a)w{z\c + a) w{x"'\h' 


-c' + 


d'- 


a)w{^"\b' -a' + d' -a) 


w{y\h + a)w{u\d + a) 




\h'- 


a)w{u"'\d' - a) 



wixiy-i.i.ux^yi 


\a' - b') w{yszi,^", ioyiz^lc' - b') 




d' — c') w^x^UijC," , XiU^l 


d' - a') 




w{x\ 


\a' + a)w{z 


\c' + a)w{x"'\b- c + 


d — a)w{z"'\b — a + d — a) 


w{y\b' + a)w{u\ 


d' + a) w{y"\ 


\b — a)w{u"\d — a) 



w{xiy3,^,ujxsyi 


\a-b) w{yszi,^",ujyizs 


\c-b) 


w{z3Ui,^,ZiU3\ 


d — c) w{x3Ui,^" ,XiU3\ 


d — a) 
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